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ABSTRACT 

Context. Star formation at earlier cosmological times takes place in an interstellar medium with low metallicity. The Large Magellanic Cloud 
(LMC) is ideally suited to study star formation in such an environment. 

Aims. The physical and chemical state of the ISM in a star forming environment can be constrained by observations of submm and FIR 
spectral lines of the main carbon carrying species, CO, C i and C n, which originate in the surface layers of molecular clouds illuminated by the 
UV radiation of the newly formed, young stars. 

Methods. We present high-angular resolution sub-millimeter observations in the N159W region in the LMC obtained with the NANTEN2 



telescope of the I2 CO J = 4 — > 3, J = 1 — > 6, and 13 CO ./ = 4 — » 3 rotational and [Ci] 3 Pi- 3 Po and 3 P2~ 3 Pi fine- structure transitions. The 
13 CO J = 4 — > 3 and [C i] 3 P 2 - 3 Pi transitions are detected for the first time in the LMC. We derive the physical and chemical properties of the 
low-metallicity molecular gas using an escape probability code and a self-consistent solution of the chemistry and thermal balance of the gas 
in the framework of a clumpy cloud PDR model. 

Results. The separate excitation analysis of the submm CO lines and the carbon fine structure lines shows that the emitting gas in the N159W 
region has temperatures of about 80 K and densities of about 10 cm -3 . The estimated C to CO abundance ratio close to unity is substantially 
higher than in dense massive star-forming regions in the Milky Way. The analysis of all observed lines together, including the [Cn] line 
intensity reported in the literature, in the context of a clumpy cloud PDR model constrains the UV intensity to about x ~ 220 and an average 
density of the clump ensemble of about 10 5 cm 4 , thus confirming the presence of high density material in the LMC N159W region. 
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1. Introduction 

Since the formation of the first generation of stars in our 
Universe, the cyclic process of formation and destruction of 
stars has progressively enriched the interstellar medium (ISM) 
with heavy elements and dust. This implies that, at earlier cos- 
mological times, stars formed in lower metallicity, lower dust- 
to-gas ratio molecular gas. Newly formed massive stars illu- 
minate their progenitor molecular clouds with far-ultraviolet 
photons (FUV; 6eV < hv < 13.6eV) producing ph oton- 
dominated regions (PDRs; iHollenbach & Tielensl 1 19991 and 
references therein). The lower dust-to-gas ratio reduces the 
ability of the low-metallicity molecular cloud to attenuate FUV 



photons, making low-metallicity PDRs more extended than 
their Solar metallicity counterparts. This has an impact on 
the C + /C/CO transition lay er in PDRs, as CO is mor e ef- 



ficiently photo-dissociated (Ivan Dishoeck & Bla ck 1988), en- 
hancing the C + and C abundance relative to CO. Therefore, 
star formation in low metallicity gas differs from that in so- 
lar metallicity gas, as the structure and thermal balance of the 
progenitor molecular gas is different. The extent to which the 



retically studied (e.t 


Malonev & Black 1988; Leaueux et al. 


1994 


: Bolattoetal. 


19991 Kaufman et alJ 199? Rollig et al. 



20061) and t o some extent obse rvationally constrained (e.g. 
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Israel 1997 ; Madden et al. 2006). The study of the properties 
of the low metallicity molecular gas is of great importance in 
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order to understand the evolution of youngest galaxies in our 
Universe. 

Nearby dwarf and irregular galaxies are examples where 
we can observe massive star formation taking place in low- 
metallicity molecular gas. One of the best studied examples 
is the L arge M agellanic Cloud (LMC). Due to its proximity 
(50kpc:lFeasdf l999) a nd ne arly face-on orientation (i = 35°; 



120-35Cfl(in units o f the intensity of the mean interstell ar radi- 
ation field derived by |Prainell978T) . lBoratto et al.l d2000t) found 



van der Marel & Ci oni 2001), the LMC is ideal to study the 



impact of m assive star format ion in its low-metallicity (Z ~ 
0.3-0.5 ZM: IWesterlundlll997l) molecular ga s. The LMC has 
been entirely mapped in [Cn] 1 58 jum bv iMochizuki et aD 
Jl994h and in 12 CO J = 1 -> bvlCohen et alJ dl988l) at low 



resolution of 8.8 arcmin. iFukui et all (fl999) revealed the first 
spatially resolved image of the giant molecular clouds over the 
LMC in the 12 CO / = 1 — » emission at 2.6 arcmin resolu- 
tion (~ 40 pc) by using the NANTEN 4m telescope. The subse- 
quent CO studies with NANTEN offer the details of the molec- 
ular and star-formation properties of these GMCs in the LMC 
dMizuno et alj|20qik lYamaguchi et aT]|200U Imitz et all 12007 : 
Fukuill2002l 120071). The CO observations of the LMC show 



that the molecular gas towards the molecular clou d complex 
associated with the N159 Hn region dHenizdl 19561) is an out- 
standing star forming region which is part of the most signifi- 
cant molecular concentration in the galaxy. An extensive study 
of the CO emission with improved angular resolution has been 
carried out within the ESO-SEST Key P rogramme on CO in 
the Magellanic Clouds dlsrael et alJI 19931 and subsequent pub- 
lications). 

Observations of low-/ 12 CO and 13 CO transitions with 
improved angular resolution (~10pc) have unveiled three gi- 
ant molecular clouds ( GMCs); N159E, N 159S, and N159W 
d Johansson et all [19981: iBolatto et~al]|2000l) . At similar scales, 
the TC nl distribution in the N159 complex has been mapped by 
iBoreiko & Betzl dl99ll) and llsrael etaD dl996l) . At lower angu- 
lar resolution of 1 10", and hence averaging over larger spatial 
scales of about 50 pc, the [Ci] 3 Pi- 3 P and 12 CO / = 4 -> 3 
transitions have been detected in this region wi th the AST/RO 
telescope dStark et al.ll 1 997t IBolatto et alj|2000b . 



The brightest CO peak in the entire LM C is N159W. In 



this re gion, an excitation analysis made by iJohansson et al 
dl998l) . based in low-/ 12 CO and 13 CO transitions, suggested 
gas densities of n — 10 4 cm~ 3 and = 10 - 20 K. A sim- 



ilar analysis by iMinamidani et al.l (120071) . including the 12 CO 
/ = 3 — > 2 transition, resulted in similar gas densities but 
somewhat higher temperatur es (/kin > 30 K). I ncluding 12 CO 
/ = 4 — > 3 observations, IBolatto et alJ d2005l) argue that the 
CO intensities are best explained by a warm (7\; n = 100 K) 
and low-density (n =10 2 crrT 3 ) and a cold (7\in = 20 K) 
and dense (n - 1 5 cm " 3 ) c omponents. Note tha t the work of 



Johansson et alJ d!998l) and IBolatto et alJ d2005l) apply to dif- 



ferent spatial scales (lOpc and 26 pc, respectively). 



The PDRs in N159W h ave been studied by Israel et all 



19961) . IBolatto et alJ d2000l) . and iPak et al.l (119981) . Based on 
different t racers of th e stren gth of the incident FUV radia- 
tion field, Israel et al. ( 19961) derived a value between x - 



that the /[ciij/^co ratio is enhanced compared to that in sim- 
ilar Galactic star-forming regions while the /[cij/^co remains 
simil ar. This is i n agre ement with a simple analytical model 
bv IBolatto et al. (1999) that predicts growth of the [C n] and 
[C i] emitting layers with decreasing metallicity, resulting from 
a shrinking of the CO molecular core. The N159W reg ion was 
included in PDR model calculations by lPak et al.l dl998l) . which 
constrained the physical conditions of the line-emitting gas to 
10 3 - 7 cirT 3 < n < 10 4 - 7 cirT 3 and 60 < x < 600. 

In this paper, we present high-resolution (lOpc) observa- 
tions of the N159W region obtained with the new NANTEN2 
telescope. We present detections of the 12 CO / = 4 — > 3, 
J = 7 -> 6, and 13 CO / = 4 3 rotational and [Ci] 3 Pi- 3 P 
and 3 P2- 3 Pi fine-structure transitions. The 13 CO J — 4 — > 3 
and [C i] 3 P2- 3 Pi transitions are detected for the first time in 
the LMC. The 12 CO / = 7 -> 6 line has been detected at the 
low angular res olution of the AST/RO telescope towards the 
30 Dor region dKimll2006l) . The NANTEN2 observations are 
described in Section 2. We derive physical properties of the 
low-metallicity molecular gas using an escape probability ra- 
diative transfer code independently for the CO and [C i] lines 
(Section 3), and a self-consistent solution of the chemistry and 
thermal balance of the gas is given using a photon-dominated 
region (PDR) model for a clumpy interstellar medium (Section 
4), explaining the observed emission of all observed lines, in- 
cluding low-/ CO lines and [C n] 158 /mi from the literature. 
We discuss our results in Section 5, and summarize in Section 
6. 



2. Observations 

We use the new NANTEN2 4m telescope situated at 4865 m 
altitude at Pampa la Bola in northern Chile to observe the 12 CO 
/ = 4 -> 3 (461.0408 GHz) , / = 7 6 (806.6517 GHz), 
and 13 CO / = 4 -> 3 (440.7654GHz) rotational and [Ci] 
3 Pi- 3 P (492.1607GHz) and 3 P 2 - 3 Pi (809.3446GHz) fine- 
structure transitions toward N159W. The line parameters de- 
rived from the NANTEN2 observations of LMC N159W are 
listed in TableQ] 

The observations were made toward the peak 12 CO / 
I -> intensity position in N159W, located at a = 5 h 39 m 36 s .8 
and 5 = -69°45'31'/9 (J2000), using the total-power observing 
mode. We use a reference position at a — 5 h 38 m 53^2 and 6 = 
-69°46'28'.'8 (J2000), which is free of 12 CO / = 1 -> emis- 
sion (3 cr upper limit of 6.2 K kms -1 ). The duration of each 
beam switch cycle was varied between 20 and 30 seconds, de- 
pending on the atmospheric stability during the observations. 
The final spectra result from total integration times of 20 mins 
to about 1 hour. The pointing was checked regularly on Jupiter, 
IRC+10216, and IRc2 in OrionA. The applied corrections were 
always smaller than 20", and usually less than 10". 



Ilsrael et all ( 119961) o riginally deriv ed this value in unit s of the 
iHabind dl968h field. The lDrainel Jl978l) field is related to the lHabind 
(1968) field by a factor of 1.71 when averaged over the full FUV 
range. 
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The observations were conducted with a dual-channel 
460/8 10 GHz receiver installed for verifying the submillimiter 
performance of the telescope. Double-sideband (DSB) receiver 
temperatures were ~ 250 K in the lower channel and ~ 750 K in 
the upper one. The intermediate frequencies (IF) are 4 GHz and 
1 .5 GHz, respectively. The latter IF allows simultaneous obser- 
vations of the 12 CO ./ = 7 — > 6 line in the lower and of the 
[Ci] 3 P2- 3 Pi line in the upper sideband. These two lines are 
observed simultaneously with one of the lines in the 460 GHz 
channel. As backends we used two acousto-optical spectrome- 
ters (AOS) with a bandwidth of 1 GHz and a channel resolution 
of 0.37km s 1 at460GHz and 0.21 km s 1 at 806GHz. 

The LMC observations with NANTEN2 were performed in 
several sessions during the September to November 2006 pe- 
riod. In all sessions both channels of the dual frequency test- 
receiver were active, so that the determination of the atmo- 
spheric transmission for both channels can be done in a sin- 
gle step with a common atmospheric model. This new atmo- 
spheric calibration scheme, implemented as part of the standard 
observing software kosma_control developed at Universitat zu 
Koln and in s talled at NANTEN2, uses a parametrization of the 
Pardo et al.l Q2004J) atmospheric model ATM. The single, free 
parameter, fitted to the observed atmospheric emission spec- 
trum derived from hot/cold/sky-measurements in the standard 
calibration cycle, is the precipitable water vapor pwv. The pa- 
rameterized model provides frequency lookup tables of the dry 
(a v ) and wet (b Y ) contribution to the opacity of the atmosphere 
above the observing site, t v = a v + b v pwv. The a v and b v 
-coefficients are derived from the ATM model atmosphere cal- 
culated for the 4865 m altitude of the NANTEN2 observa- 
tory. The cold-load is a liquid nitrogen load, the hot-load is 
at the cabin temperature, which is monitored electronically. 
The sky temperature is taken as the ambient outside temper- 
ature monitored by the weather station of the observatory. This 
assumption is occasionally verified through sky-dip measure- 
ments, which are also used to determine the forward/spillover 
efficiency. Using this calibration, the relative intensities of lines 
observed in different receiver bands simultaneously is, beyond 
the signal/noise of the data, limited only by the precision and 
validity of the atmospheric model. The main advantage of this 
new scheme is that it gives a much higher precision for the 
determination of the atmospheric line-of-sight transmission in 
receiver bands with low transmission and poorer sensitivity, be- 
cause it includes the stronger (hot-sky) difference signal of the 
higher transmission bands and the higher signal-to-noise data 
from the more sensitive bands in determining the transmission 
at each frequency point. It also fully takes into account fre- 
quency dependent variations of the atmospheric opacity across 
the receiver bands, such as resulting from narrow atmospheric 
lines and tails of wide ones. Comparison with the traditional 
calibration approach which determines the atmospheric trans- 
mission as an average over a representative section of each 
receiver band individually from hot/cold/sky-measurements, 
shows a much higher robustness against errors in the values 
for the ambient and effective atmospheric temperatures used in 
the calibration. 

The other factors that enter in the relative calibration of the 
intensity of lines from different receiver bands are the ratio of 



their main beam efficiency and the pickup from the error beam 
pattern outside the main beam. Main beam effic iencies were 



determined from Jupiter radio continuum scans dSimon etal 



2007), assuming a brightness temperature of Jupi ter of 154 K 



and 1 26 K respectively at 460-490 and 810 GHz (Griffin et al 



1986). Comparison of the fitted angular width of the Jupiter 



scan with the disk diameter of Jupiter allows to derive the 
FWHM of the main beam, and from its brightness, the main 
beam efficiency. The thus derived observational parameters are 
listed in Table [2] A full description of this new calibration 
scheme for ground based submm observations is in preparation 
for a separate, technical paper. 

A coarse estimate of the error beam of the NANTEN2 tele- 
scope was obtained by analyzing scans across the edge of the 



ing 
al. 



Moon and the Sun (iSimon et alJ 120071) . These show a signal 
consistent with an error beam size of some 400" at 460 GHz 
and 240" at 810 GHz, consistent with the diffraction pattern of 
individual panels, which have a size of 50 by 70 cm, and hence 
due to panel mis-alignment. We estimate the integrated power 
(relative to the main beam) of these error beams to be 10% and 
15% at 460-490 GHz, respectively 810 GHz. The error beam 
intensities are consistent with a surface rms due to panel mis- 
alignment of about 25 /mi, as independently confirmed by pho- 
togrammetric measurements. From the analysis of these data, 
we estimate the ratio of main beam efficiencies to be deter- 
mined to a precision of about 10%, resulting in a similar sys- 
tematic error in the observed intensity ratios of lines in the 460 
and 810 GHz band. 

The simultaneous coverage of both the 12 CO ./ = 7 — > 6 
and [C i] 3 P2 — > 3 Pi transitions in the upper and lower sideband 
implies that their relative calibration, including the atmospheric 
imbalance between the two sidebands, is very precise, limited 
only by the radiometric S/N in this receiver channel. The lower 
frequency channel of the dual channel receiver was alternat- 
ingly tuned to the 12 CO J = 4 -> 3, the [C i] 3 Pi - 3 P transition, 
and the 13 CO J = 4 — > 3 line. The atmospheric transmission 
was particularly good for an initial set of observations covering 
the 12 CO J = 4 — > 3 line in the lower channel. At this time the 
zenith opacity was about 0.40 at the 12 CO / = 4 -> 3 line fre- 
quency and about 0.92 for the 12 CO J = 7 — > 6 line, resulting in 
very high S/N spectra in particular for the 12 CO J — 4 — > 3 line . 
Using the a v and b v coefficients from the IPardo et al. (2004) 



model these zenith optical depths correspond to a precipitable 
water vapor (pwv) of about 0.5 1 mm. A subsequent observing 
session covering the [Ci] 3 Pi- 3 Po transition had again a very 
good zenith opacity of about 0.5 to 0.6 at 492 GHz, and about 
0.7 to 0.8 in the upper frequency channel for the CO 7 — > 6 
line (pwv = 0.38 to 0.44 mm). During the 13 CO 7 = 4^3 
line observations the weather was very good for a short ses- 
sion, with an opacity of about 0.5 to 0.6 at the line frequency of 
440 GHz, at the edge of the 600 pm atmospheric window, and 
correspondingly low zenith opacity of about 0.6 at the 12 CO 
J — 7 — > 6 frequency (pwv =* 0.31 mm). Most of the 13 CO 
./ = 4 — > 3 data were taken in a long integration at substantially 
higher zenith opacity, varying around 1.2-1.5 (corresponding to 
about 1.4 to 1.5, pwv 0.70 to 0.89 mm, at 806 GHz). While 
this implies that the absolute calibration at these high opacities 
may easily be influenced by systematic errors (such as non- 
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Fig. 1. LMC N159W: observed lines and Gaussian fits (red 
line). 

perfect temperature estimates for the effective radiative tem- 
perature of the atmosphere), the relative calibration between 
the low- and high-frequency line is largely unaffected and thus 
rather accurate. 

In the following, we thus used estimates of the systematic 
errors for the line intensity ratios introduced through the cali- 
bration uncertainty (as explained, these are mainly due to the 
precision with which the beam efficiency ratio is determined), 
as quoted in Tableland add these quadratically with the radio- 
metric noise, assumed to be given by the formal errors of the 
Gaussian fits. 

3. Observational Results 

3.1. NANTEN2 data 

The spectra observed towards LMC N159W are shown in 
Fig. [1] together with Gaussian fits of the center, width and 



integrated line intensity, resp. amplitude. The results of the 
Gaussian fits are listed in Table Q] The line center and width 
determined from the high signal-to-noise 12 CO J = 4 — » 3 
spectrum were taken and held fixed for the Gaussian fit of the 
line amplitude of the other lines. This is adequate as their pro- 
files are, within their limited signal-to-noise, consistent with 
that of the I2 CO 7 = 4^3 line. 

As explained above, due to the simultaneous observations 
and the joined atmospheric transmission calibration for the 
lines in the upper and lower frequency channel, the intensity 
ratio of the lines observed is rather precisely determined and 
is mainly limited by the estimated systematic errors due to 
the determination of the beam efficiencies as discussed above. 
The only further correction concerns the different beam size 
of the low and high frequency channel of 38", resp. 26". 
Unfortunately, we did not have enough nights with superb 
weather to make small maps around the peak positions at (0,0) 
in order to spatially smooth the high frequency observations 
to the angular resolution of the low frequency beam. For cor- 
recting the beam size effects for the high frequency observa- 
tions we thus have to rely on an independent determination of 
the source extent. AST/RO has repeatedly observed the 12 CO 
J = 4 — > 3 line, once w ith a relatively large beam FWHM of 
204" (Bolatto_etaL|[2000), resulting from improper alignment 
of the AST/RO optical components, giving an integrated in- 
tensity of 8.4 K kms -1 , and with the diffr action beam size o f 
the AST/RO telescope at 460 GHz of 109" dBolatto et alJ2005h . 
giving an integrated line intensity of 20.8 K km s _1 . Assuming a 
Gaussian source distribution with FWHM © s and a source peak 
brightness temperature of T S)Pea k, the beam filling correction 
gives a main beam brightness temperature in a beam of FWHM 
&i, of T„,t, = @ 2 /(® 2 + ©?)r, iPeffl t . The ratio of the main beam 
brightness temperature in two different beams 1 and 2 is given 
by Rx, 2 = T^i/Tnbi = (®% 2 + ®?)/(®^ + ©?)• With this, we 
can derive the source size from the observed brightness ratio in 
two beam sizes: 



© 



#1,2©; 



6,1 



•®L 



1 -#1,2 



giving ® s = 52" from comparing the NANTEN2 12 CO J = 
4 — > 3 data with the early, and 0, = 61" for comparison with 
the later AST/RO data. 

We thus adopt a source size of 60" and correct the observed 
intensity ratios for the beam coupling of the different beam 
sizes at two observing bands x and y to the source intrinsic 
intensity ratio R s via 



R, = 



T mb ,y 0?+0 2 / 



This correction assumes that the source structure is the 
same at the spatial scales observed with the 460 and 810 GHz 
channels. The validity of this assumption likely dominates the 
uncertainties in the line ratios. However, as we will see in 
Section|4]and|5] the good fit obtained for the physical parame- 
ters supports this assumption. Table[3]shows the line ratios thus 
corrected and used for the excitation analysis discussed below. 
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Table 1. Gaussian-fit line parameters derived from the NANTEN2 observations of LMC N159W and from the literature 



Line 


dllipillUUC 


L-CllLCi 


r vv jnivi 


i n tf^ m"Q 1 
lllLCgl al 




^K^ 


[JsJll a J 


[Kill a J 


[ l\. Kill a J 


NANTEN2 data 










12 CO J = 4 -> 3 


9.31 (0.05) 


238.3 (0.016) 


7.32 (0.04) 


72.6 (0.3) 


12 CO J = 7 -> 6 


3.33 (0.09) 








13 C0 7 = 4-> 3 


0.95 (0.10) 








[Cl] 3 P, -> 3 P 


1.51 (0.05) 








[Ci] 3 P 2 -^ 3 P, 


1.87 (0.08) 








literature data w 










12 CO J = 1 -> 


6.79 (0.15) 


238.0 (0.073) 


6.95 (0.17) 


50.0(1.6) 


13 coy = 1 -> 


0.78 (0.08) 


238.0 (0.35) 


6.81 (0.83) 


5.6 (0.9) 


[Cn] 


7.7 (0.9) 


232.8 (0.8) 


9.8 (0.9) 


75.5 (7.0) 



<n) amplitude derived from fitted Gaussian width and integral 
{b) center and width fixed to 12 CO J = 4 — » 3 best fit values 



w references for literature data see Section 3.2 
Table 2. NANTEN2/dual-channel test receiver: Observational Parameters 



receiver channel 


460-490 GHz 


810 GHz 




0.5 


0.4 


beam FWHM ["] 


38 


26 


Tree [K] 


250 


750 



3.2. Additional Lines from the Literature 



Beyond the lines newly observed with NANTEN2 we have 
included in the analysis described in the following also the 
12 CO and 13 CO J = 1 -» lines (the 12 CO data are pre- 



sented in Pin eda et al.1 d2008l) and the 13 CO data were kindly 
provided by Jiirgen Ott), as w ell as the [C nl 1 5 8 pm line ob- 
served spectrally resolved by iBoreiko & Beta (1199 II) to bet- 
ter constrain the physical parameters of the emitting gas. The 
12 CO and I3 CO / = 1 — * lines have been observed with the 
ATNF Mopra 22 m telescope fully-sampled with its 33" beam 
at 115 GHz and later convolved to 45" resolution, and the in- 
tensities are converted to the main beam brightness temperature 
scale using rj m \, = 0.42. The [Cn] 158 //m line has been ob- 
served with the Kuiper Airborne Observatory with its 43" beam 
( Boreiko & Betzfl 99 1 ) at 158 /mi and its intensity is quoted on 
a main beam brightness temperature equivalent intensity scale. 
The pointing of the KAO observations is about 25" further 
south than the NANTEN2 position. We take their values with- 
out further corrections. In fact, the integrated intensity of the 
[C n] line from the KAO observations of 4.3 x 1CT 4 erg crrr 2 s _1 
sr~ 



4. Escape probability radiative transfer: excitation 
analysis 

We now use the observed line intensity ratios and absolute in- 
tensities to constrain the physical properties of the emitting gas. 
We start by using only the intensity ratios, which are indepen- 
dent of the beam filling of the source. We analyze separately 
the CO emission and the atomic carbon emission, but assume 
that both originate from the same material in the beam with 
identical beam filling factors. This is motivated by the fact that 
the beam sizes at 460-490 and 810 GHz at the distance of the 
LMC corresponds to linear scales of 10 and 7 pc, respectively. 
Moreover, both CO and [C i] emission arise from dense molec- 
ular material, more specifically from the UV illuminated sur- 
face of cloud structures unresolved at this linear scale. 

We use the escape proba bility radiative transfer model by 



Stutzki & Winnewisser dl985h to calculate the line intensities 
as a function of kinetic temperature cloud density tijn and 
column density per velocity interval N/Av of the species of 
interest (note that the column density parameter in the radia- 
tive transfer model is the radial column density of the spheri- 
cal clump; the beam averaged column density in a beam larger 



is consistent within 20% with the integrated intensity read tna n the clump thus is a factor of \ {^f) / (nr 2 ) _ f higher) 



off the [C n] map published by Israel et al. ( 1996b . observed ve 
locity unresolved with the FIFI spectral imaging instruments on 
the KAO. 

To estimate adequate errors for the line intensity ratios, we 
assume that the absolute calibration of all the independently 
observed lines is about 20%. Table [3] quotes the thus derived 
total errors, obtained by quadratically summing up the formal 
fit error and the calibration uncertainty. 



The collision rate coefficient s used are from Flower & Launav 
(1 19851) for CO-H 2 and from lSchroderet all dl99l|) for C°-"hT 
The escape probability model produces peak brightness tem- 
peratures of an individual spherical model clump, averaged 
over the projected area of the clump on the sky. To compare 
with the observable brightness of the model clump, this has 
to be scaled down by the beam filling factor of the individual 
clump. Within the standard model of a molecular cloud consis- 
tent of many fragments which move around with a virialized 



6 Pineda, Mizuno, Stutzki, Cubick et al.: Submillimeter Line Emission from LMC N159W 

Table 3. Observed line integrated intensity ratios, corrected for beam coupling to an assumed common source size of FWHM of 
60" 



Species 


Ratio (fl) 


fit error 


calibration error 


total error 


12 CO J = 7 -> 6 / CO J = 4 -» 3 


0.30 


0.008 


10% 


0.031 


[Ci] 3 P 2 -^/[Cl]^ -* 3 P 


1.05 


0.06 


10% 


0.12 


13 CO 7 = 4^3/ 12 CO J = 4 -> 3 


0.10 


0.011 


10% 


0.015 


[Ci] 3 P 2 -* 3 Pi /CO J = 7^6 


0.56 


0.052 




0.052 


12 CO / = 1 -> / 12 CO J = 4 -> 3 


0.69 


0.10 


20% 


0.17 


n COJ = 1 ->0/ l2 COJ = 1 -»0 


0.12 


0.02 


20% 


0.03 


[Cn]/ 12 CO J = 4-+3 


1.1 


0.1 


20% 


0.24 



velocity dispersion, the ratio of the observed to the model line 
intensity, once corrected by the beam filling factor, is thus given 
by the ratio between the clump internal linewidth and virial ve- 
locity dispersion. For optically thin emission of the individual 
clumps, and also for moderately optically thick clumps, the ob- 
served/model intensity ratio gives the ratio of beam averaged 
column density to clump intrinsic column density, as long as 
the total integrated line intensity still increases linearly with 
the number of clumps (no velocity crowding). We will use this 
below, when we compare the absolute line intensities of 12 CO, 
13 CO and [Ci] with the model predicted intensities. 

4.1. [Ci]-emission 

The observed line intensity ratio of the two [C i] lines of about 
1 .05 (Table [3]) is compared with the escape probability radia- 
tive transfer prediction in Fig. [2] We present model line ratios, 
in a stack of 3 plots corresponding to kinetic temperatures of 
40, 80, and 120 K, as a function of density and radial column 
density per velocity interval, covering the range from 10 2 to 



10 6 



cm in density, and 10 to 10 cm / (kms ) in col- 



umn density per velocity interval! The line ratios and absolute 
intensities discussed further down are presented in the same 
way and cover the same parameter range. Clearly, the observed 
value can only be reproduced by the model at moderately high 
temperatures of above about 60 K, and at densities of the [C i]- 
emitting material of at least 2 x 10 3 crrT 3 at high temperatures, 
and slightly higher, around 6 x 10 3 cirT 3 at the lower limit of the 
temperature range. At these densities, slightly higher than the 
critical densities of the atomic carbon fine structure lines, the 
carbon lines should be close to thermalized. Fig. [2] also shows 
the model predicted absolute intensity of the [Ci] 3 Pi — > 3 Po 
line. The analysis of the different CO-lines presented below ar- 
gues for a beam filling factor of about 1/6, so that the observed 
peak brightness of about 1 .5 K should be scaled up to an intrin- 
sic brightness of about 9 K corresponding to a C i column den- 
sity model parameter of 10 17 cm _2 /(kms _1 ). In comparison to 
the kinetic temperature range derived above and together with 
the densities inferred, we conclude that the low observed and 
beam filling corrected brightness implies that the [C i] lines are 
well in the optically thin regime. Thus, multiplication by the in- 
ferred beam filling factor is appropriate. The intersection of the 



line ratio and absolute line intensity curves, lying on the verti- 
cal branch of the [C i] line ratio curves in the density/column 
density plane, independent of the value of the beam filling fac- 
tor, implies that the density of the emitting material is in fact 
close to the lower limit values quoted above. 

4.2. CO-emission 

^.Frorn the NANTEN2 data we ob- 
tain two CO line intensity ratios, namely 
r mb ( I2 CO/ = 7^6)/7 mb ( 12 CO./ = 4^3) and the iso- 
topomeric ratio r mb ( 13 CO / = 4 -> 3)/r mb ( 12 CO / = 4 -> 3). 
The predicted line intensity ratios are plotted in Fig. [3] This 
plot (and the following ones) assume a fractional abundance 
of 13 CO/ 12 CO of 40 by shifting the contours for the ' 3 CO 
lines accordingly along the column density axis. A value for 
the I3 CO/ I2 CO abundance ra t io in the N 159W regio n of 40 
is cited by iJohansson et al. | dl994h and llsrael et alj J2003) 
quote a value of 50, whereas iHeikkila et al J ( 19991) derive a 
substantially lower value of only 21. Their abundance ratios 
for other 12 C and 13 C carriers, namely H 13 CO + and 13 CS, 
are higher. This may indicate, that the isotopomeric CO ratio 
is indeed affected by fractionation, whereas the elemental 
12 C/ 13 C ratio is around 40. Fractionation may occur due to the 



exchange reaction 13 C + + 12 CO 



3 CO + 35 K being 



2 Note that the beam average column density per velocity interval 
of an unresolved clump, as explained above, is a factor of 4/3 larger. 



exothermal and thus preferentially incorporating 13 C into the 
molecular form. Selective photodissociation of 13 CO due to 
the lower self-shielding would cause the opposite. PDR model 
calculations show that the proportion of these effects strongly 
depend on the temperature derived for the C + /C i/CO transition 
zone and the comparison between different models gives 
non-conclusive answers. If we constrain the density to values 
compatible with those derived from the [C i]-analysis above, 
the two observed CO line ratios can only be reproduced by 
the model at a temperature of about 80 K. The assumption of 
equal densities is fully justified, because their is no physically 
plausible scenario that has CO at low densities of about 100 
cm -3 , while the Ci is located at higher densities of close to 
10 4 cm 4 , as long as we assume that the lines originate in a 
PDR-like structure. Note that the reverse, i.e. [C i] originating 
from much lower density material than CO, namely from a 
more diffuse halo of dense molecular clouds, would be also 
reasonable. 

The intersection of the two observed CO-line ratios con- 
strains the CO-column density model parameter to a value of 
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about 10 17 cm~ 2 /kms _1 at 80 K. Recall that we show the plots 
for an abundance of 13 CO/ 12 CO of 40. The column density 
thus derived changes roughly in proportion to the assumed 
isotopomeric abundance ratio, as is confirmed by running the 
model for different values; we do not show the corresponding 
figures here, in order to save space. 

Fig. |4] shows that the model predicted clump averaged ab- 
solute line intensity of the 13 CO/ = 4 — > 3 line at the column 
densities and in the density regime constrained by the observed 
line ratios as discussed above, is about 6 times higher than the 
observed one at a temperatures of about 80 K. This implies a 
beam filling factor of the CO material of about 1/6. 

We now consider in addition the line intensity ratios in- 
volving the 12 CO and 13 CO J — 1 — > transition. The line 
intensity ratio between 12 CO / = 1 -> and 4 -> 3 of 0.69 
traverses the density/column density regime of n = 10 4 cm -3 
and about 10 17 cm 2 / (kms 1 ) derived from the other CO line 
ratios only at temperatures of above around 80 K and higher 
(see Fig. [5]). The isotopic J = 1 — > line ratio does not traverse 
this regime at all. It always requires substantially higher col- 
umn densities than allowed by the line ratios of the higher-/ 
lines. This implies that reducing the 12 CO/ 13 CO abundance ra- 
tio from the value of 40 assumed here to lower values, will 
bring the / = 1 — > ratio into better agreement, although Fig. [5] 
shows that lowering it by just a factor of 2 down to about 20 
would not be sufficient. We report (without showing a figure), 
that the agreement becomes reasonably good at lower tempera- 
tures of 20 to 40 K, indicating that a temperature gradient with 
the lower-/ CO emission originating in lower temperature gas 
and showing stronger isotopic fractionation can give a consis- 
tent picture, th us confirming the low 12 CO/ 13 CO ratio of about 



20 reported by Heik kila et all dl999l) . 



To summarize, the escape probability radiative transfer 
analysis implies that the CO and CO emission observed in 
N159W, if it originates in gas with densities close to 10 4 cm" 3 
as constrained by the [C i] emission, comes from gas with ki- 
netic temperatures of about 80 K and with a source intrinsic 
total column density of 10 17 cm _2 /(kms _1 ) and a filling factor 
of about 1/6. The 12 CO / = 7 -» 6, / = 4 -> 3, / = 1 -> 
and the 13 CO / = 4 — » 3 and / = 1 — > line ratios taken 
together are not consistent with a single temperature compo- 
nent, but require a range of temperatures and possibly isotopic 
fractionation. 

4.3. CO and [Ci] column densities 

The escape probability radiative transfer analysis gives column 
densities for both CO and Ci of 10 17 cm~ 2 /kms _1 . But the 
12 CO and 13 CO / = 1 — > lines indicate an additional, cooler 
component. The derived column density ratio [Ci] / CO of 1.0 
thus is an upper limit due to the unknown contribution of the 
lower temperature CO-material. 

As discussed above, the beam average column density is 
a factor of 4/3 higher than the column density model parame- 
ter. Together with the observed line width of 7.3 kms _I and 
the beam filling factor of 1/6, both CO and Ci thus have a 
beam average column density of 1.6 x 10 17 cm~ 2 . Assuming 



the canonical abundance CO/Fb of 8 x 10 5 and the same for 
C i, both tracers correspond to a mass in the beam of 3000 M 
each, and 6000 M in total. Considering that our beamsize of 
38" FWHM corresponds to 9.5 pc at the distance of the LMC 
and hence covers only part of the N159W cloud's extent of 
about 15 pc, this is in good agree ment with the CO-lumi nosity 
mass of 1.9 x 10 4 M derived bv LFohansson et all dl998l) . 



5. PDR-model analysis 

In the previous section we have analyzed the excitation con- 
ditions of both [Ci] and CO/ 13 CO separately. These condi- 
tions correspond to typical values of the line-emitting mate- 
rial. In the following, we investigate whether the observations 
can be explained in terms of PDRs distributed in a clumpy in- 
terstellar medium. The modest temperatures, with a tendency 
towards lower temperatures for the CO lines than for the [C i] 
lines, are roughly consistent with those expected for the mid-/ 
CO and [C i] -emitting layers in a PDR illuminated by the 
modest UV- field of x about 175, independently derived for 
the LMC N159W region from the FIR continuum luminosity 
dlsrael et al. 1996). Compared to typical star forming regions 
in the Milky Way, the estimated column density ratio of C i rel- 
ative to CO of up to unity is rather high and this is expected 
to be the case in a lower metallicity environment such as the 
LMC with z=0.4. Similarly, the [Cn] emission is expected to 
be stronger due to a larger extent of the C-ionized PDR sur- 
face layer resulting from the weaker UV attenuation by dust 
an d weaker self shieldin g of CO at low metallicities. 

Bolatto et al. dl999h . in a simple semi-analytical model 



(their model A) predict this behavior for the carbon species 
with the assumption of finite size clumps making up the emit- 
ting gas; the PDR depth scale is determined by the UV absorp- 
tion per nucleon, which is due to heavy elements (in the form 
of dust, carbon etc.) and thus, referred to the hydrogen column 
density, scales up inversely with the metallicity. The [C n] and 
[C i] emitting layers thus grow in size inversely with metallicity 
and the molecular core of the clump, where the CO emission 
originates, shrinks. 

The KOSMA-t PDR model dStorzer et all 1 1996b allows 
to self-consistently calculate both the chemical and temper- 
ature structure of spherical clumps and con firms these ba- 
sic trends with metallicity (Rolli g et al. I l2006h . but takes into 
account all details such as the variation of the temperature 
in the emitting region of each species with varying metal- 
licity. In order to model the emission from a cloud com- 
plex such as the N159W region, one additionally has to 
take into account the fact, that the emission is unlikely to 
arise in a single clump, but that the emitting complex shows 
structure on all scales. Quantitative studies of the struc- 
tural pro perties in many star forming regions in the Milky 
Wav (e.g. IWilliams et alii 1994 lElmegreen & Falgarone|l996 : 
I Stutzki et al.lll988HKramer et all 1 9981: iHeithausen et all 1998 : 
Hever & Terebevl 1 19981) . but also on larger scales in the 
Milky Way and external gala x ies dStanimirovic et aL 2000; 



l200lt 



Kim et al 



Stani mirovic & Lazarianl 1200 It iFukui et al 
120071) . show that the observed velocity channel maps can be 
decomposed into clumps with a power law mass spectrum 
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dN/dM oc M~ a with a clump mass spectral index of 1.8 to 
1.9; without the velocity dispersion due to turbulence, this 
would not be possible, as in many cases the clumps overlap 
along the line-of-sight and can only be separated by the ad- 
ditional velocity information. An independent approach is to 
stud y the power spectra of the ob served maps, either directly 
(e.g. lErmegreen & Falgaronel 1996b or via the A-variance anal- 
ysis dStutzki et all 1998b. which shows that these often follow a 
power law. IStutzki et alJ (119981) showed, that both approaches 
are consistent with the additional assumption of a power-law 
mass-size relation for the individual clumps, M oc r y , with an 
index y of about 2.3. Note that values of y > 2 imply increas- 
ing density with decreasing clump size, i.e. denser small scale 
structures, as expected for a turbulent medium, where the small 
scale structures are created by local compression due to tran- 
sient waves; y = 2 would imply constant column density of 
the structures, independent of size. Although the demonstra- 
tion of a clear mass-size power law and the derivation of the 
exact power exponent is rather difficult and uncertain, in the 
few cases where it has been possible it confirms values around 
the one quot ed above; an example is the analysis of the Polaris 
flare region dHeithausen et al.|[l998 ). 

Based on this scenario, an obvious approach is to model the 
observed emission of a cloud complex by combining the emis- 
sion of many clumps, individually modeled with a clump PDR- 
model such as KOSMA-r, which follow a power law mass 
spectrum and mass-size relation. This approach has been suc- 
cessfully applied to several i ndividual regions dKramer et al.1 
l2008t iMookeriea et al.ll2006l iKramer et alj|2004 l2005h . but 
also to the global emission of the Milk y Way in the submm /FIR 
cooling lines, as observed by COBE ( Cubick et al.l l2007). We 
will use it here to analyze the combined set of observed lines 
in LMC N159W. T he details of the model are explained in 



Cubick et al. (2007). 



Fixing the clump mass spectral index and the mass-size 
power-law index to a = 1.8 and y = 2.3 respectively, the 
clumpy cloud complex is characterized by the mean density 
of the total clump volume (each clump being modeled as a 
Bonnor-Ebert sphere with oc r~ 1,5 density law, so that the clump 
average density is twice the surface density), the total mass of 
the clump ensemble, and the upper and lower mass cut-off of 
the mass spectrum. As long as the individual clumps are not 
too densely packed spatially and in velocity, so that their emis- 
sion does not overlap along the line-of-sight and in frequency, 
changing the total clump mass simply changes the absolute 
brightness of the clump ensemble by adding or removing ma- 
terial with identical mass distribution; equivalently, it changes 
the area and volume filling factor of the clump ensemble. In the 
model fit, it is thus used as the parameter that adjusts to the ab- 
solute intensity observed and it can be ignored, as long as one 
considers only relative intensities. 

Keeping the other parameters fixed, but changing the up- 
per mass cut-off shifts the distribution of the material to larger 
or smaller clumps and thus will affect the relative intensities 
of the various lines according to their preferential origin from 
the clump cores or surfaces. Lowering the lower mass cut-off 
simply adds more and more smaller clumps, which, as long 
as the lower mass cut-off is sufficiently low, are small enough 



to correspond basically to PDR-surface layers; given the above 
value of the mass-spectrum index, the small clumps do not con- 
tribute to the total mass, i.e. amount of material, and thus this 
change does not significantly affect the observe d emission. A 
somew hat open issue, discussed in more detail in lCubick et al.l 
d2007l) concerns the fact, that one can argue about the stabil- 
ity and life-time of such small clumps and to what extent the 
mid-7 CO emission from these small and dense structures af- 
fects the overall emission of the ensemble in these lines. We 
do not discuss this issue further here. For the models discussed 
below, the low mass cut-off was fixed to 10~ 3 M . 

The free model parameters are thus, as long as we consider 
only relative intensities, the UV-field, parametrized by x m 
units of the Draine field, the clump-volume mean density n, the 
metallicity z and the upper mass cut-off M max . For the follow- 
ing analysis, we first fix the metallicity to the canonical value 
for the LMC of z = 0.4; later on we discuss, how the results 
vary with metallicity. We then do a least-square fit in the n 
-parameter plane for different values of M max . The least-square 
fit analysis includes all line ratios listed in Table [3] i.e. the rel- 
ative intensities of the full suite of low-/ and mid-7 CO and 
I3 CO lines, both [Ci] lines, and the [Cn] line. Unfortunately, 
the PDR model grid is presently only available for a 12 C/ 13 C 
elemental abundance of 60, slightly higher than the experimen- 
tally determined value for the LMC N159W cloud. 

We get the best fit, with a normalized xt in = 6 for a UV- 
field of^- around 220 and a density n of a few times 10 5 crrT 3 . 
This is obtained at a rather low upper-mass cut-off of 0.1 M . 
Fig|6] shows the ^-distribution in the parameter plane and the 
contours of the observed line ratios overlayed. Considering that 
we fit seven line ratios with 2 parameters, this is an excellent 
fit. 

We next discuss the dependency of the best fit param- 
eters on varying the upper mass cutoff and the metallicity. 
Representatives-contour plots are shown in Figj7] Increasing 
the upper-mass cut-off to 1 M and 10 M increases lhexf ni „ to 
a value of 7.5 and 12.8. This shallow rise of less than a factor of 
2 inx 2 per decade in M max , shows that the upper mass-cutoff is 
not very stringently constrained, but indicates that the clumpy 
cloud PDR scenario requests the emitting gas to show lots of 
surface and to be broken up into small scale structures. Varying 
the metallicity to the solar value of 1, the^, is 15 and 22 re- 
spectively for a M max of 0.1 and 10 M , so that the PDR model 
clearly favors the canonical metallicity of the LMC in order to 
reproduce the observed intensity ratios. Fig|7] shows that the 
discussed variations do not much alter the best fit value for the 
UV-intensity and the density. 

Considering the uncertainties about the carbon isotope el- 
emental abundances to use for the LMC and the open issue 
of isotope selective fractionation in PDRs, it is worthwhile to 
consider a PDR-model fit excluding the I3 C-carrying species. 
Without showing the corresponding plot, we report that the 
best-fit density and UV-intensity are unaffected; this is due to 
the fact, that the line ratios involving the 13 C-carrying species 
vary rather shallow across the parameter region of interest. A 
close inspection of Fig|6]in fact shows, that neglecting the ob- 
served [C n]/CO J = 4 —> 3 ratio in determining the x 2 - 
minimum would allow a solution at slightly lower densities 
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of about 10 5 crrr 3 and higher UV- intensities of slightly above 
^-=1000, where all other lines intersect. 

Whereas the value for the UV-intensity is consistent with 
that independently derived from the FIR luminosity, as dis- 
cussed in the introduction, and thus only confirms this result 
independently, the density regime favored by the clumpy-cloud 
PDR-model, n — 10 5 cirT 3 or slightly above, is clearly much 
higher than other estimates of the density in the LMC N159W 
region. Together with the small value for the upper-mass cut- 
off this argues for a very clumpy structure with high density 
clumps to explain the observed line ratios. 

In the last step of the analysis, we can now adjust the pre- 
dicted absolute line intensities to the observed value, by chang- 
ing the total mass of the ensemble to fit the absolute inten- 
sity of e.g. the 12 CO J — 4 — > 3 line. This gives a mass in- 
cluded in the beam of 2.8 x 10 4 Mp, in good agreemen t with the 
CO-luminosity mass quoted bv lJohansson et al.1 dl998l) and the 
mass derived above from the escape probability analysis. Note 
that the clumpy-PDR model fits all lines simultaneously, where 
as the low-/ CO line analysis is biased towards the lower ex- 
citation material and the escape probability radiative transfer 
analysis above emphasized the warm and dense gas traced by 
the mid-7 CO transitions and the [C i] fine structure lines. 

We conclude that the emission of all observed CO lines, 
from J = 1 -> to 7 -> 6 for 12 CO, J = 1 -> and 4 -> 3 for 
13 CO and both [C i] and [C n] fine structure lines, is consistently 
modeled as originating from a dense, highly clumped medium 
illuminated by an FUV-freld of x * 220 in the LMC N159W 
complex. The model fit constrains the densities to values of 
about 10 5 cirT 3 and preferences the canonical metallicity of 
the LMC of about z = 0.4. 



6. Discussion 

Previous CO excitation analyses suggest a low-temperature 
(T - 10 - 20 K) molecu lar gas in N159W (iJohansson et al 



compared with what is typically found in Galactic dense 
massive st ar formation region s (0.04-0.2; for a compila- 



1998; iBolatto et al.1 120051) . However, these studies are only 



based on CO observations which trace more shielded material 
in the PDR. The excitation analysis using the [C i] fine structure 
line intensity ratio allow us to better explore the conditions at 
the cloud surfaces, as th is ratio is a very good tracer of the sur- 
face temperature (e.g.lStutzki et alll997l) . IBolatto et all d2005l) 
and Minamidani et al. (2007) derived larger temperatures for 
this region (T = 150 K and T > 30 K, respectively). Our work 
confirms a high density of the bulk of the CO-traced molec- 
ular material in the N159W region. Wide ranges of possible 
gas densities have bee n sugg ested by [Johansson et al.l (|l998) 
and iMinamidani et al.1 d2007h (TO 3 - 10 s cm" 3 ). Although 
poorly constrained, IBolatto et al. I (f2005h also proposed a two- 
component solution of low density (« = 10 2 cirT 3 ), high tem- 
perature (T = 100 K) and high density (n = 10 5 cm~ 3 ), 
low temperature (T = 20 K). High densities of the more 
shielded material are also suggested by observations of molec- 
ular species li ke CS, SO, and H 2 CO (n = 10 5 - 10 6 cirT 3 ; 
Heikkilaetal1ll999l) . 

We derive a [Ci]/CO column density ratio toward N159W 
close to 1.0, possibly slightly lower, depending on the con- 
tribution from a cooler CO component. This value is large 



Mookerjea et al. 



tion, see 

larger C abundance relative 



2006), and therefore argues for a 
to CO as a result of the en- 



hanced CO photodissociation in low-metallicity PDRs. Note 
that the [Ci l/CO column density ratio is not constant in 
the Galaxy dMookeriea et al. 20061) . High values are also 
found in the less dense outsk irts of massive star-forming 
regions (Mookeriea et al. 20061) and in high-latitude clouds 
(Inga lls et all 1 19971) . Recent observations have also shown 
high [Cil/CO column density ratios in extragalactic nuclei 
(Hit schfeld et al.ll2007l) . In order to determine the relative im- 
portance of the main carbon species, we derive the C + col- 
umn density, assuming Mem = 6.4 x 10 20 /[cirj crrr 2 (erg s~' 
cirT 2 sr dCrawfordetalJ [l985rT and us ing the /rem = 
4.3 x 10~ 4 erg s~' cirT 2 sr 1 observed by Boreiko & Betzl 
dl99ll) . With this method, we obtain Afcm = 2.75 x 10 17 



crrr . Using this value and the CO and CI column densities 
of /Vco.ci = 1-6 x 10 17 cm -2 as derived above, we obtain a 
relative abundance of the main carbon species in N159W of 
Afcrrj : Afcrj : /V C o=46%:27%:27%. In Galactic star-forming 
regions CO is the most abundant carbon species. For example, 
in the DR21(OH) star-forming region, the relative abundance 
of the main carbon species is Mcni^Mc n : 'Vco=3%:10%:87% 



(Jakob et al. 2007). As we can see, the abundance of C + and 



C relative to CO are indeed enhanced in the low-metallicity 
gas in N159W. The different distribution of carbon is reflected 
in the enhanced /[cii]// 12 coj=i^o an d Aci]/^ 12 coy=i^o li ne ratios 
observed in LMC N159W, relative to Solar metallicity PDRs. 
^From the data used in this paper, we obtain /[cii]//i 2 coj=i^o 
1.13 and /[cipPi- 3 p // 12 coy=i^o - 0.22. These values can be 
compared with the ratios deter mined from COBE all-sky mea- 
surements ( Wright et al Jl 1 99 ll) of 0.97 and 0.13, respectively. 

Polycyclic aromatic hydrocarbons (PAHs) are known to 
have an important influence in the carbon chemistry, increas- 
ing th e C and decreasing the CO abundances in clo ud sur- 
faces dLepp & Dalgarnolll988l; iBakes & Tielenslll998b . In or- 
der to assess any effect of PAHs in the PDR modeling, we 
also compare the observations with a single-clump version 
of the KOSM A-t model that include s PAH in the chemical 
network (see iPineda & Benschl d2007l) for details on the in- 
clusion of PAHs in the KOSMA-t model). P AHs were also 



consid ered in the PDR model calculations by iKaufman et al 



dl999l) . We find that the /[crj//co line ratios typically increase 
by up to a factor of 10 when PAHs are considered. As a re- 
sult, in order to match the observations, even larger densities 
than derived above would be required (n = 10 6 - 10 7 cirT 3 ). 
The effects of PAHs in the chemical network of the clumpy- 
PDR model will be explored in a future work. The PAH abun- 
dance might be reduced in low-metallicity molecular gas, as 
observations revealed that P AH emission is weak or absent in 
low metallicity galaxies (e.g.lMaddenl20 00; Houc k et al. 2004; 
Eng elbracht et al.ll2005l;lMadden et al.l2006l) . However, the ob- 
served PAH emission spectra in the N 1 59-5 Hn region rese m- 
bles the spectra of Galactic H n regions dVermeii et al ] |2002|) . 



3 This is an approximation for an excitation temperature of > 91 K 
and a density larger than n a = 5 x 10 3 , cm" 3 
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The strength of the FUV radiation field required by the 
PDR model to explain the observations (x — 220 ) is consis 



tent w ith the independent determination made by llsrael et al 



1996), using UV, radio continuum, and IRAS 60yum emis- 



sio n (y = 120 - 35 0). Using the same assumptions made 
by llsrael et al. (1996) to determine x from the IRAS 60 /im 
data, but using the high-angular resolution Spitzer 70/im map 
( Mei xner et al. 2006), we obtain a consistent value ofx = 177. 
This strength of the FUV radiation field might be consid- 
ered low for a massive star form ation region such as N1 59W 
(for example, in the Or ion Bar, Wal mslev et al. | (120001) and 
Young Owl etai] d2000l) find that x - 10000). However, 
Israel et al. (1996) argue that the low metallicity and low dust- 
to-gas ratio of this region increases the FUV photon free path 
lengths, resulting in a greater geometric dilution of the FUV 
radiation field. 

7. Summary and conclusions 

In this paper, we derived the physical conditions of the line- 
emitting gas in the N159W region in the LMC. We compared 
an escape-probability code and a PDR model with NANTEN2 
observations of the 12 CO J = 4 -> 3, J = 7 -> 6, and 13 CO 
7 = 4^3 rotational and [Ci] 3 Pi- 3 P and 3 P 2 - 3 Pi fine- 
structure transitions. Our results can be summarized as follows: 

- The analysis of the excitation conditions for both the CO 
submm-lines and the [Ci] fine structure lines shows tem- 
peratures of about T — 80 K and densities of about n — 10 4 
cm -3 for the emitting gas. 

- The estimation of the C and CO column densities suggest 
a high C/CO abundance ratio of close to unity compared to 
that found in Galactic dense massive star-forming regions. 
This is an indication that CO is more readily destroyed 
in the low-metallicity environment of N159W. Comparing 
with an estimation of the C + column density, we confirm 
that the distribution of the main carbon species in PDRs 
is altered in the low-metallicity molecular gas, with larger 
C + and C abundances relative to CO. This result is in 
agreement with observations of an enhanced /[cin/^co and 
7[d]//co ratios in N159W, compared with values observed 
in solar metallicity PDRs. 

- The PDR modeling, calculating the emission as resulting 
from an ensemble of spherical clumps with a power law 
mass spectrum and a power law mass-size relation, provid- 
ing a simple representation of the fractal structure of the 
ISM, gives a consistent fit to all observed line ratios for a 
moderate strength of the FUV radiation field of ~220^o- 
This value is co nsistent with a pre vious independent esti- 
mation made by llsrael et al. ( 1996b . It constrains the av- 
erage density of the clump ensemble to about 10 s cm -3 , 
confirming previous indications of high-density CO-traced 
molecular material in the LMC N159W region. The clumpy 
cloud PDR model also gives a relatively small upper clump 
mass cutoff for the clump ensemble, indicating a very much 
broken up structure of the ISM in N159W. 

Similar studies in other regions in the LMC with 
NANTEN2 combined with observations from upcoming obser- 



vatories such as Herschel and SOFIA, will improve our knowl- 
edge about the impact of star formation in the low-metallicity 
molecular gas. This will help us to improve our current un- 
derstanding of the processes of star formation in the young 
Universe. 
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Fig. 3. Comparison between observed and escape probability 
Fig. 2. Comparison between observed and escape probability model predicted 12 CO 7 = 7— >6to/ = 4— > 3 (solid contours) 
model predicted [C i] line ratios (dashed contours) and absolute and 13 CO to 12 CO J = 4 — » 3 (dashed contours) line ratios. The 
intensity of the [Ci] 3 Pi — > 3 Po line (solid contours). Contour plots assume a fractional abundance of 13 CO/ 12 CO of 40. 
levels are given at the top in the format "min(step)max", with 
the minimum contour level plotted as a thicker line. The levels 
of the dotted contours are the observed values. For the abso- 
lute intensity contours these are shown additionally scaled by a 
beam filling factor of 1/3 to 1/6. 
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Fig. 4. Escape probability model predictions and the observed 
absolute intensity of 13 CO J = 4 — > 3. The three dotted 
contours show the observed value and the observed intensities 
scaled up by a factor of 3 and 6 respectively, to match the col- 
umn density derived from the line ratios at the temperature of 
about 40 K, respectively 80 K. 



Fig. 5. comparison between observed and escape probability 
model predicted 12 CO J = 1 —> and J = 4 -» 3 line ra- 
tios (dashed contour) and the 13 CO to 12 CO / = 1 -> line 
ratio (solid contour). 
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Fig. 6. Clumpy KOSMA-r PDR-model fit to the observed 
line ratios overlaid onto the gray shaded contours of the cor- 
responding reduced^ 2 . Note that the observed line ratios are 
plotted as solid lines, with the lower and upper limits as dot- 
ted respectively dashed lines. The x 2 gray-scale ranges from 
one to four times the minimum x 2 -value of 6 at a density of 
log(« [cm -3 ]) = 5.2 and an incident FUV-flux of \ogix) =2.1. 
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Fig. 7. Reduced x 2 contours of the clumpy KOSMA-r PDR- 
model results for different values of upper clump mass limit 
Mcimax of the clump ensemble and metallicity Z. Note that all 
contours are given in values of the individual minimum x 2 at 
the corresponding values of the density and FUV-flux inserted 
in the legend. 



